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Cotton flow tends to evolve a given initial metric on a three manifold to a conformally flat 
one. Here we expound upon the earlier work on Cotton flow and study the linearized version 
of it around a generic initial metric by employing a modified form of the DeTurck trick. 
We show that the flow around the flat space, as a critical point, reduces to an anisotropic 
generalization of linearized KdV equation with complex dispersion relations one of which is an 
unstable mode, rendering the flat space unstable under small perturbations. We also show 
that Einstein spaces and some conformally flat non-Einstein spaces are linearly unstable. 
We refine the gradient flow formalism and compute the second variation of the entropy and 
show that generic critical points are extended Cotton solitons. We study some properties of 
these solutions and find a Topologically Massive soliton that is built from Cotton and Ricci 
solitons. In the Lorentzian signature, we also show that the pp-wave metrics are both Cotton 
and Ricci solitons. 


I. INTRODUCTION 

Gravity has a smoothing effect on a mass or matter distribution, as is evident from all the 
sufficiently massive spherical objects in the universe, including our planet which has deviations 
from surface smoothness with its tallest mountain and deepest point in the ocean amounting 
to less than 10“^ of its radius. For more massive compact objects, the deviations from surface 
smoothness would be even less: For example on a neutron star with 1.4 times the mass of the 
Sun, the height of a mountain could be around 1 cm (not much to hike!). Black holes are even 
smoother: They have "no hair". Rotations about an axis can introduce larger deviations from 
perfect spherical shape, but that is a totally different story. For the earth, this deviation is about 
4 X 10“^, that is the polar radius is around 25 km smaller than the equatorial one. While gravity’s 
smoothing effect on objects that we can see is clear, what we also know is that space-time, even 
if devoid of any matter, as a manifold itself, also gravitates due to the non-linearity of the field 
equations. This raises many questions that fall in the intersection of physics and mathematics. For 
example, one might wonder if one can study the topological classification of manifolds with the 
help of gravitational equations of some sort. The answer is affirmative but not obvious and actual 
execution of such a classification is certainly not straightforward. First, one needs to make sure 
that there is a connection between the Riemannian structure (the metric, the curvature) and the 
topology. After all the metric and the curvature are local quantities but global topology a priori 
should not be determined by the local structure. This cursory look is actually not correct and there 
is a connection between the topology of a manifold and its Riemannian structure as beautifully 
exemplified by the formula relating the Fuler number of a closed two surface in terms of its total 
Gaussian curvature: 27rx(F) = d?x. 

This and analogous results in higher dimensions encourage one to study the global topology of 
manifolds with the help of curvature and the metric. Gonsidering the above mentioned smoothing 
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nature of gravity, and the fact that topological classification boils down to identifying manifolds 
that can be smoothly deformed to each other, one might perhaps try to let the metric on a manifold 
evolve in time with Einstein’s equation or some modified form of it. While this seems like a feasible 
idea, it is also well-known that singularities, such as black hole singularity, form in a sufficiently 
general non-linear theory of gravity. Hence one must somehow deal with singularities or better yet 
turn them into advantage as is done in Ricci flow with surgery. 

The choice of gravity-like equations in order to study the topology of manifolds is an important 
issue. For example, if one naively takes the Einstein’s equations (say in vacuum), then depending 
on the number of dimensions, one has a very different nature of solutions. Considering three 
manifolds per se, one can split the four dimensional Einstein’s equations into 3 -|- 1 as space and 
time hence study the evolution of 3 spaces in time. But it turns out that there are all sorts 
of solutions, including propagating local modes (gravitons), gravitational waves etc. due to the 
hyperbolic nature of the equations. Sticking to three dimensional manifolds, the goal is to have a 
finite dimensional solution space (up to perhaps scaling of the metric and the diffeomorphisms). In 
the choice of the proper equations, uniformization theorem of two dimensional closed, orientable 
surfaces, which reduces to finding conformally equivalent metrics of constant Gaussian curvature 
(usually normalized as 1,0,—1 for the sphere, torus and higher genus tori respectively), gives us 
an important hint: That is, one must find "constant curvature" metrics on a given manifold. Here 
by "constant-curvature", we mean either constant scalar curvature or constant Ricci curvature or 
in the extreme case, constant Riemann curvature. Of course the last one is too restrictive and 
its solutions are already known; they are the maximally-symmetric spaces. On the other hand, 
finding constant Ricci curvature metrics boils down to solving Einstein’s equations in vacuum 
{Rij = ^Qij) for which a general strategy is not known even though there are books compiling 
plethora of solutions. The question whether a manifold of dimension larger than two admits a 
constant scalar curvature is called the Yamabe problem [l|, which in some sense generalizes the 
two dimensional uniformization problem to higher dimensions. 

All the above discussion suggests that given a differentiable manifold M and an initial metric 
go on it, we should try to get a "nicer", more symmetric metric by smoothly deforming the initial 
metric. The way to proceed is clear: As in the case of the usual diffusive scalar heat equation, one 
must have an equation of the form 


dt 


gij{t,x) = Eij{t,x), 


( 1 ) 


where the deformation parameter t is "time" (or a parameter related to the energy scale) and it is 
an external parameter not a dimension on the manifold, while x denotes the local coordinates on 
the manifold which we shall take to have a positive signature metric. Depending on the property 
of the manifold one is interested in, a specific Eij tensor can be chosen. For example, to define 
the Ricci flow, which was introduced by Hamilton Q], one takes the following "volume-normalized 
" equation 

(2) 


where r is the average Ricci scalar 


r = 


(^-x^/g R, 


where V = 


d^x^. 


( 3 ) 


Fixed points of the flow are Einstein metrics, therefore this is in some sense a parabolic extension 
of Einstein’s equation, albeit its parabolic nature is not apparent at this non-linear level. Namely 
with this flow, the metric is deformed in such a diffusive way that the manifold eventually admits a 
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constant Ricci curvature while preserving volume. Hamilton used this equation to give another 
proof of Poincare’s uniformization theorem for two dimensional closed surfaces mentioned above. 
Over a decade ago, Perelman 0,0, using a modified form of the Ricci flow equatiom proved 
the three dimensional uniformization theorem, then called the Thurston’s conjecture [y] with a 
heroic effort bringing Hamilton’s expectation and programme to a successful end. Even though 
the uniformization theorem itself and the equations employed to prove it are so closely related to 
physics, not much work has been done on Ricci flow in physics (save a couple which we mention 
below). 

Let us expound upon how the uniformization theorem and Ricci flow type equations are related 
to physics. First of all, just about the same time (and actually, a little before) Hamilton introduced 
the Ricci flow, the equation was found by Friedan 0 as the renormalization group (RG) equations 
of the couplings in a non-linear two dimensional sigma model. In the non-linear sigma model, the 
metric is just like any other coupling and runs with energy, hence there is a /3 function for it and 
the (5 function is given as the Ricci tensor plus infinitely many two tensors built from the curvature 
tensor and its derivatives. Two-dimensional sigma model fields are coordinates on the target space 
which is a Riemannian manifold. RG flow in the sigma model corresponds to a modified Ricci flow 
in the target space. 

One can also see the relevance of the three dimensional uniformization theorem in general 
relativity (GR): In the 3-1-1 dimensional formulation of GR, three manifolds evolve in time. Since 
closed 3 manifolds are classified by the uniformization theorem, in essence, the theorem is at the 
heart of the question of what the shape of the apparent 3 dimensional spatial universe is. Let us 
give another example of how these flows can be relevant to physics. Consider a gravity theory (a 
Euclidean one) defined with certain boundary conditions. There could be many solutions satisfying 
the boundary conditions. In that case, the question is which solution is the global minimum of the 
Euclidean action. One can use a Ricci flow type equation to study the stability of the solutions and 
determine possible transitions between solutions. In fact, this was done in 0 for general relativity 
in a box with a boundary of x with 3 critical points in four dimensions. Their exercise is a 
demonstration of how Euclidean version of gravity (which could possibly correspond to quantum 
processes such as black hole decay in the Lorentzian setting) can be studied by geometric flows. 
Namely, the discrete solution space is extended by the Ricci flow hence an off-shell theory can be 
constructed. See 0 for a similar analysis using the Ricci-Cotton flow that studies the transition 
between the two vacua Anti de Sitter (AdS) and warped AdS. See [l0(], for some other possible 
applications of Ricci flow and 1]| for a modified Ricci flow with additional fields and for some 
higher derivative flows. 

Gotton tensor {Cijj 


Inspired by the work on Ricci flow, in 
flow exclusive to three dimensions as 


1 ^ 


was used to define a geometric 


dt 


gij{t,x) = KCij{t,x), 


(4) 


where fixed points are locally conformally flat metrics. Note that since the Weyl tensor vanishes 
identically in three dimensions, a manifold is conformally flat if and only if Cij = 0. Therefore, 
Cotton flow equation could be used to answer the question: whether a given three manifold Ai 
admits a conformally flat metric or not. In this sense pure Cotton flow is orthogonal to the 
Yamabe flow since Yamabe flow keeps the conformal class intact, while the Cotton flow changes 
the conformal class. In 0 , equation ([3]) was derived as a gradient flow of an entropy functional 
and the equation was numerically and analytically solved for homogeneous Thurston’s geometries. 
In particular, it was shown that a deformed initial metric on flows into the conformally flat 
round metric on S^. In the table below, we compile these homogeneous geometries and compute 
their Cotton tensors just for the sake of completeness. As it is clear from the table, the first five 
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of the geometries are fixed points of the Cotton flow. Moreover, homogeneously deformed initial 
metrics on 5^ x and x R^ are left intact under the Cotton flow. On the other hand, 

under the Ricci flow 3“^ x R^ and x R^ geometries are degenerated 1^. Analytic solutions 
show that homogeneous metrics in the Nil class tend to a pancake degeneracy and the Solv metrics 
develop cigar degeneracies. 


Thurston Geometries 


Geometry 

Metric 

Curvature Tensors 

Cotton Tensors 


= dx'^ - 1 - dy'^ + dz^ 

Rij = 0 

Cij = 0 


ds^ = dx"^ + sin^ xdy"^ -\- {dz -I- cos x dy)'^ 


Cy = 0 


ds^ = ^{dx‘^ + dy‘^ -I- dz^) 

Rlj = 

Cy = 0 

52 x Ri 

ds^ = dx'^ -I- sin^ x dy^ -|- dz^ 

Rii = 1 , R22 = sin^ x 

Cy = 0 

X Ri 

ds^ = ^{dx‘^ + dy‘^) -t dz'^ 

Rii = R22 = —^ 

Cy = 0 

Sol 

ds^ = e^^dx'^ + e^'^^dy^ -t dz^ 

R33 = -2 

C12 = 2 

Nil 

ds^ = dx'^ - 1 - dy"^ + [dz -I- x dy)"^ 

D _ Id _ 1 

-rtii — 2? -^22 — 2 

R33 = 5, R23 = 1 

Cll = -\,C 22 = X^-\ 
C33 = 1 , C23 = X 

SL{ 2 ,R) 

ds^ = ^{dx^ + dy^) + {dz + ^dy)^ 

p 2 _ 19 p _ 3 

^ij 45-^ 2’ -^11 2^’ 

R22 = -yj, R23 = R33 = ^ 

Cii = C22 = 

C23 = - 1 , C33 = -2 


In this current work, we shall study several pertinent points in Cotton flow, the most important 
being the linear (in)stability of flat space. We also show the instability of Einstein spaces and 
some conformally flat spaces. We shall also linearize the Cotton flow equations about a generic 
background and apply DeTurck’s technique [1^ to get a slightly better equation, yet we shall 
not be able to prove short-time existence for an arbitrary initial metric. Nevertheless, ensuing 
discussions along this line will need our linearized equations. We also refine the earlier gradient 
flow formulation of the Cotton flow and define new objects which we shall call extended Cotton 
solitons that are critical points of the entropy functional. We study two examples of Cotton and 
Ricci solitons and define Topologically Massive solitons. There has been several other interesting 
works in the literature reg arding the Cotton flow (or its more complicated cousin Ricci-Cotton 
flow). For example, see 16l4l8l| for Cotton-like flow in the Horava-Lifshitz gravity context. See also 
H-E] where Cotton solitons were found. 

The lay-out of the paper is as follows; In Sec.II, we linearize the Cotton flow equation about 
a generic background and identify the symmetries of the equation. In Sec.Ill, we fix the diffeo- 
morphisms and scalings of the metric via the DeTurck trick to simplify the equation and compute 
the principal symbol of the relevant operator. In Sec.IV, we study the linearized stability of the 
flat space which is a critical point of the flow and find an unstable mode. In Sec.V, we refine 
the gradient flow formulation of the Cotton flow, and expand the entropy functional for a generic 
fixed point of the flow up to second order in perturbation theory, and show that Einstein spaces as 
critical points are unstable. In Sec VI, we give details of the properties of extended Cotton solitons 
and show that there are no non-trivial compact solutions. We give examples of Cotton and Ricci 
solitons as well as Topologically Massive solitons. In the Appendices A,B,C, we expound upon 
some of our calculations and in Appendix D, we show how the Cotton tensor flow under arbitarary 
flows and specifically the Ricci flow. 
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II. LINEARIZATION OF THE COTTON FLOW ABOUT AN ARBITRARY CURVED 

BACKGROUND 


Cotton flow on a 3 manifold is defined as 


1 


dt9ij — 


(5) 


where we have scaled t to set k = 1 and from now on we suppress the arguments of the tensors. 
The Cotton tensor is given as 


Akl 




( 6 ) 


with being the anti-symmetric tensor density defined as = -|-1. [In Appendix A, we compute 
how some other tensors and scalars built out of tensors flow under Cotton flow.] In what follows 
we will use the anti-symmetric tensor 




ijk 


V9 


(7) 


is symmetric and covariantly conserved = 0, and traceless gijC^^ = 0. These constraints 

tell us that out of 6 components of the (symmetric) metric gij, only 2 can be determined by the 
Cotton flow equation ([S]), hence a gauge fixing must be done. To understand the local nature of 
the Cotton flow about a flat background, equation Q was linearized for flat metrics gij = 5ij in 
0. In this section, we will generalize this result for a generic background gij to understand the 
linearized version of the flow. The results are relevant to the existence and uniqueness of the flow 
as well as the linear stability of the critical points, namely whether the critical points are saddle 
points, minima or maxima of the action from which Cotton flow is derived. For the existence and 
uniqueness of the flow, unfortunately we have nothing to say since the equation is a third order 
PDF and necessary mathematical technology does not seem to exist yet, but stability issue will be 
studied in the next section. 

To linearize the Cotton tensor about an arbitrary background gij, it is useful to recast it in an 
explicitly symmetric form as 2^ 


2cb = 


( 8 ) 


where G^i = Wi — iR is the Einstein’s tensor. Assuming gij = gij + hij, and hij is small 
compared to the background metric gij which at this stage arbitrary (namely not a critical point 
of the flow), then one obtains the linearized Cotton tensor as 


2(C*^')l = - y R^kh^i 

111 1 0 

- S^iVkh - h^iVnS^k + Snk^^h^i + g'^' Si^^nhh + i^j, 


(9) 


where Sij = Rij — \gijR is the traceless Ricci tensor and □ = g^^'^i^j and h = g^^hij. Note that 
"L" refers to the linearization and all the barred quantities are taken with respect to the background 
metric, which also raises and lowers the indices of linearized quantities. One must be careful with 
the linearization of the up and down indices, for example 


{Gij)L = {G^^)L9ik9ji + G'^jhik + G^hjk. 


( 10 ) 


^ One is almost tempted to call this "100 % Cotton flow," as there will be non-pure, mixed flows that we shall discuss. 
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Derivation of Q is somewhat tedious, hence we relegate the details into Appendix B. This is the 
main equation that should be used in the study of the stability of various critical points of the flow. 
Let us note several properties of Q. First of all one obtains the contraction gij{C'^^)L = —hijC^^, 
as expected from the linearized tracelessness condition {gijC'^^)L = 0. Secondly, the linearized 
version of the divergence-free condition of the Cotton tensor reads at this order as 

Vi{C^^)L + {V\k)LC'^^ + {TUk)L&^ = 0 , ( 11 ) 

which is satisfied by Q which can be shown after a rather long computation. 

Let us also understand the symmetries of the linearized Cotton tensor Q. Under infinitesimal 
diffeomorphisms and scalings of the metric as 

= ViCj + VjCi + ^{x)gij, (12) 

with A(x) a scalar function and Q a vector field, Q transforms as 

(13) 

where the first three terms come from the Lie derivative £^(C'*-^)l as expected. The important 
point is that (lliip vanishes at the fixed points (C*-^ = 0) of the flow, hence diffeomorphism and 
scalings apparently make the critical point a saddle point. But this is actually a red-herring coming 
from the symmetry of the theory and so must be moded out or gauged out to study genuine flows. 
This property of the linearized flow equation can also be reached in the non-linear flow equation. 
Diffeomorphism are easy to understand (they are similar to the Ricci flow case), therefore let us 
concentrate on the scalings of the metric. 

Consider changing the t-parameter to s —>■ s(t), and define a new metric on the manifold as 

gij{s) = ip{s)gij{t{s)), (14) 

such that gijit) satisfies the Cotton flow equation 

dtgijit) = Cij{t). (15) 


Then, lets check the flow satisfied by gij{s): 


d ~ r \ ~ 1 dcp(s) , xdt „ , . 

Since the Cotton tensor transforms under conformal scalings as Cij[g) 
y/gC^j is conformally invariant) one has the following flow equation 

d ^ . dClogipis)) ^ 3. .dt~ , , 

gij + ^^{s)—Cij{g). 


ds 


(16) 

V~^Cij{g), (namely 
(17) 


Hence, choosing = 1 yields 


= + ( 18 ) 

Therefore by parameterizing the time variable one can produce the scaling term. 

Next having observed the symmetry of the flow, we can introduce the DeTurck trick [l^ which 
was used to show that the Ricci flow is parabolic and hence diffusive. Without this trick, parabolic 
nature of the Ricci flow for an arbitrary initial metric is a highly cumbersome task to show [^. 
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Since the Cotton tensor is of third order, one should not expect to apply the theory of elliptic 
differential operators. But, as we shall see, a modified version of DeTurck trick removes the zero 
modes of the relevant operator at the critical point and hence modes out the symmetries of the 
flow, yielding a somewhat simpler equation. 

Before closing this section let us briefly note how the Cotton flow equation should be linearized in 
the first order formalism where instead of the metric one uses the dreibein and the spin connection. 
For this purpose, let e“ and w°‘h denote the dreibein and spin-connection 1-form respectively. Then 
the curvature 2-form reads 


= dw^b + A (19) 

from which the Ricci 1-form and the scalar curvature can be computed via the inner product. 

(Rzc)“ = R = iaibR^^"- (20) 

Then the Cotton 2-form reads 

+ w^b A y^ (21) 

with the Schouten 1-form given as Y°' = {Ric)°‘ — \Re°‘. With these definitions, the first order 
form of the Cotton flow (which was employed in [l^ ) to discuss the flow of homogeneous quantities 
reads 


5 ^ 6 “ = ( 22 ) 

where * denotes the Hodge dual. Assuming e“ satisfies 

= 0, (23) 

following the notation of [^, let us expand around the critical point of the flow as 

e“ = e“ + (24) 

Here 7be“ = 5^. The linearized flow equation reads 

e\dtip\) = (25) 

To compute the linearized Cotton 2-form, let us first find the linearized spin connection 

w'^b = + s'" {Db^°'c - Da^Pbc) , (26) 

where Da = TaD and the background spin connection satisfying De°‘ = de“ -|- w°'b A = 0. One 
obtains the curvature 2-form as 

Rab — Rab 6 A D^Db^fac Daifbc'j: (27) 

and the Ricci 1-form as 

{Ricr = {Ricf + - DbD^ip\) + - D^^ba) - (28) 

where = DaD°^, from which the linearized Schouten and Cotton forms can be computed. 
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III. DETURCK TRICK IN THE COTTON FLOW 

By moving the □ operator to the left in the second term and playing with the relevant indices, 
we have deliberately arranged ([U]) in such a way that we can gauge away some terms by extending 
the Cotton flow as 


dtdij — Cij + Vihj + ^ jVi- (29) 

As in the Ricci flow case we can take the symmetry breaking vector as 

= (30) 

but to simplify our equation further by removing the third term and its partner in Q we will take 
it as 

R, ^ V™(rL - rim) +9^k9^^{r^n - fL), (31) 

where is a fixed connection on the manifold which we shall take it to be coming from the initial 
metric gij. There is a subtle issue here which can be seen for the flat {Rij) critical points: One 
might mistakingly think that the first piece in (18Ih would be sufficient to remove all the zero modes 
but it turns out that it only removes the diffeomorphisms, but not the scalings 5xhij = X{x)6ij. 
Linearization of (j8ip yields 

{V')l = V^Vkhm' + Vmh^^ - (32) 

Since one has dtg^^ = —(C*-^ + + V-^R*), the third term in ([9]) and its symmetric partner are 

canceled with this gauge-fixing, DeTurck-Cotton flow becomes 

dth}^ = -|- i -H- j -|- Reaction Terms, (33) 

where the reaction terms do not contribute to the principal symbol of the operator on the left-hand 
side. Of course, even after this simplification one cannot say much about the diffusive nature of 
the equation, since the principal symbol of the relevant operator vanishes as we show here. 

Let us calculate the principal symbol of the Cotton tensor as an operator on the metric using 
the above linearized form. We have 

a[C]{(p){h)ij = (34) 

with = g^^ipipj and ipi is a non-vanishing vector field. Then the inner product vanishes 

{a[C]{ip){h)ij ,hij) = 0. (35) 

This class of the operators are known as sub-elliptic operators when the operator is of second order. 
The classification for third order operators has not been done to the best of our knowledge, so we 
have not much to conclude from (I35p . 

As we noted, the critical points of the flow satisfy Cij = 0, that is, they are the locally con¬ 
formally flat metrics. Among these critical points, one clearly has Einstein spaces Sij = 0 as a 
subclass, including the flat space. In 3 dimensions, since the Riemann tensor can be written in 
terms of the Ricci tensor and scalar curvature, Ricci flat metrics are Riemann flat. In the next 
section we study the stability of flat metrics as a subclass of critical points of the flow. 
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IV. STABILITY ANALYSIS OF THE FLAT SPACE 

Let us now study the stability of the flat space in Cotton flow. This problem is interesting on its 
own since flat space is a fixed point of the flow, but it is also relevant to understand the stability of 
any manifold under small wave-length (or high momentum) perturbations. Here what we mean by 
small wave-length perturbation is that <C 1, with \Rij\ a typical curvature of the manifold in 
some coordinates. Namely for these sufficiently high energy modes, any space is locally flat. One 
expects that these modes die out, or diffuse as time passes if the critical point is a stable point 
under linear perturbations. But it will turn out that there is a potentially unstable mode in fiat 
space at the linearized level. 

For the flat background, the linear Cotton flow equation with the DeTurck terms becomes ^ 

dthij = d'^dkhji + did^hkj - ^didjh + i^ j, (36) 

where is the usual Laplacian in flat space. In a broad sense, this equation is a tensor version 
of the linearized Korteweg-de Vries (KdV) equation albeit an anisotropic one, in the sense that 
the flow of hij with given i and j is determined by hmn with m and n not necessarily equal to i 
and j. This anisotropy is extremely crucial for the diffusive or non-diffusive nature of the equation 
since the sign of hij is not necessarily positive. Therefore, in the stability analysis, the sign of the 
amplitude of the perturbation is also relevant. So to see the nature of this equation, we must first 
diagonalize it. 

Given the initial perturbation at t = 0 as hij{0,x) with the asymptotic fall off as 
lim|x|_>.oo hij{0, x) —>■ 0, we can construct “plane-wave” type solutions or Fourier modes from which 
we can construct the general solution that satisfies the asymptotic boundary condition. 

Let us insert the plane-wave ansatz into (1361) 

h,,(t,x)=C.,(p)e*(P—( pW, (37) 

where Cij(p) are (complex) polarizations which only depend on the mode (or momentum p) and 
w = rc(p) is the dispersion relation to be determined below. This insertion leads to a consistency 
condition on the polarizations 

2 i i 

iij = ~ + ^P^P^^ + i^ j, (38) 

where p = |p| and the trace of the polarization tensor is ^ We will now play with 

this equation to find w = w{p) and the constraints on the polarization tensor. [Note that these 
constraints come from the linearized Bianchi identity and the tracelessness of the Cotton tensor.] 
Taking the trace of (|38p . one obtains 

(l -(39) 

and contracting with yields 

(l —p^^py^ij = (40) 

To get the possible dispersion relations, let us assume that C / 0 and w = y then from (13^ 
P^P^^ij = 0 is obtained. But (14011 gives ^ = 0 so this is a contradiction which can be resolved by 


^ See Appendix C for another computation of the linearized instability. 
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assuming w = ip^ and ^ = 0 and p^'p^^ij = 0. This necessarily says that , that is ^ij is a 

transverse traceless tensor. Going back to one has 






(41) 


This equation can be diagonalized once the left-hand side is plugged to the in the right-hand 
side and contraction of r/’s are carried out, to get^ 





p. 




TT 

ij 


= 0 . 


(43) 


Since we cannot take to be zero (otherwise hij = 0), we have 



(44) 


which contradicts with our earlier assumption w = ip^. Hence we conclude that the assumed 
dispersion relation is not correct: w ^ ij?'. 

Clearly the same reasoning works for the other possible dispersion relations w = —2ip^ coming 
from (HHIl . Namely, two apparently possible dispersion relation in (I5UD and (HOD are do not yield 
possible solutions of (I55|) . Hence we must assume w ^ ip^ and w ^ —2ip^. Then, assuming ^ 7 ^ 0 
and p^p^Cij 7 ^ 0 equations and ([411)1 yield 


(1 + = 0, 

V w / 

which leads to another possible dispersion relation 


w = —ip 


(45) 


(46) 


and the identity ^ = p^p^^ij- Now, let us evaluate whether this mode is viable or not. For this 
purpose decompose ^ij as follows 

+ {PiPj ~ ^^ijP‘^)Mp) + (47) 

where A{p) is arbitrary at this stage. Using ^ = p^pP^ij in (14711 gives ^ = A{p)p‘^, and so one arrives 
at 

{'IS) 

Substituting (145)1 into (155)1 yields 

C + ^7 - ^p^7^ ^7 - (49) 


® Alternatively, for TT-modes, by taking one more derivative of (ll-ibl) . one obtains 

d?hij = -^idyh,„ ( 42 ) 

whose Fourier transform yields (1441) . 
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Inserting w = leads to the conclusion that for these modes, the TT polarizations vanish: 

= 0. Therefore this is a valid dispersion relation for the following modes 


iij 



(50) 


For these modes hij{t,x) = hence as time passes, these modes decay namely they are 

diffusive modes. 

Let us now consider the final possibility, that is ^ 0, then ^ = 0, p^p^Cij = 0) hence we have 
the complex cubic dispersion relation 


w = ±^P^\p\- ( 51 ) 

With this dispersion relation, we still need to make sure that there is a non-trivial solution of 
equation (1411) . We can recast that equation for both signs of the dispersion relation as a matrix 
equation 


^ = B^ + ^B, 


(52) 


with ^ as a 3 X 3 symmetric matrix and B as a 3 x 3 anti-symmetric matrix whose elements 

are given as 




(53) 


Equation (|52p is known as Sylvester’s equation which can be solved 2^ but we do not need 
the general solution, all we need is a non-trivial solution which we can find as follows. Suppose 
p = (0,0,p), then, ^ and the B matrix read as 


62 0 
^ = I 62 -61 0 
0 0 0 


B = 



Equation yields 


62 = 


(54) 


Normalizing .^n = 1, ^12 = depending on the sign of the dispersion w = ±|p^|p|. Therefore 
we have non-trivial solutions for both signs of the dispersion relation. 

Since the equation (j36p is linear, then the most general solution reads as 


hij{t,x) = j cfp Cij{p)e' 

We can recast this in terms of the initial perturbation as 

hij{t,x.) = J d^xG{x,x,t) hij{0,x), 

where 


ip.x±ip2|pp 




,ip.(x-x )±ip2|p|t 


(55) 


(56) 


(57) 
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This integral can be evaluated for the minus sign, in terms of hyper-geometric and Bessel-kind 
functions. On the other hand, the positive sign mode is problematic. It is clear that the negative 
sign modes decay in time and hence these perturbations are diffused. As for the positive sign modes, 
there is a linear instability. But bear in mind that this analysis has been at the perturbative level 
and the perturbation theory breaks down for these growing modes hence a strict conclusion cannot 
be made with regard to the non-linear instability of flat space in Cotton flow. 

In fact, even though we have not been able to prove stability at a non-linear level, we conjecture 
that these growing modes could be tamed at the non-linear level. If this turns out not to be the 
case, then, since the above analysis is also valid for any manifold for high momentum fluctuations 
that see the manifold locally flat, any critical point of the flow would be unstable, leading to the 
conclusion that the theory has no minima but just saddle points. This would be highly unexpected. 


V. COTTON FLOW AS A GRADIENT FLOW AND THE ENTROPY FUNCTIONAL 

In 0 a gradient flow formulation of Cotton flow was given which we refine it here to better 
understand the role of symmetries and the Cotton solitons and also study the issue of stability 
using the entropy functional. The Chern-Simons action on a manifold without a boundary 

+ ?r”r&). (58) 

has the first variation under an arbitrary change of the metric as 

5J^ = [ (59) 

Jm 

Observe that for 6gij = ViVj + VjVi + X{x)gij, the first variation vanishes, 6J^ = 0, owing to the 
fact that Cotton tensor is traceless and covariantly conserved. Hence such a variation gives us the 
symmetries of the theory: which are diffeomorphisms and arbitrary scalings of the metric. To get 
a steepest descend, we can choose the following flow 

^Qij = dtgij = Cij + ViVj -h VjVi -h X{x)gij, (60) 

which leads to 

bT=^=l d^x^gCifi^K (61) 

J M 

To get a gradient flow we must require ^ > 0 where the bound is saturated for the critical points 
that satisfy the following equation 

Cij -|- VjV)' -|- VjV) -|- \{x)gij = 0, (62) 

which we shall call it as the extended Cotton solitons. [In 0 , the case of constant A metrics 
were called Cotton solitons]. It is clear that (I6ip vanishes at the critical points, since by definition 
^9ij = 0 there. Requiring the volume density to be a conserved quantity of the flow dt^/g = 0, 
yields, from (I5U|) . the constraint 


2ViU* + 3A(x) =0, 


(63) 


for all metrics gij that are solutions of the flow. 
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To understand the stability of the fixed points, let us expand the Entropy functional of the 
Cotton flow up to the second order in 6gij for a fixed time t. Then from (15911 . one obtains 

+ {C^^)LS9iX (64) 

Jm ^ 

where gij refers to the critical points that satisfy (I62|) . Making use of ([U]) one arrives at 

d^x^g!^ - hi.VkXphh + hi.V^VkVnhr 

- hSi^Vkh^j - hi.Sk^V^Ki (65) 

+ hi.SnN^hK + hi, SkrN^h^i + , 

where we have dropped some boundary terms and made use of the Lichnerowicz operator action 
on a symmetric two tensor as 

Xi^hij = —Ohij — 2Rikjih^^ + ( 66 ) 

which yields in three dimensions 

hij = —D/iij — 2gij Rjfih^^ — 2hRij + gijhR — Rhij + SRi^hj^ + SRjkhi^. (67) 

To ensure the stability of the critical points one must have 




> 0 . 


9 ij 


( 68 ) 


For a generic conformally flat fixed point, or for Cotton soliton backgrounds, even though it is 
possible that 5‘^J- is positive, we have not been able to show this. Instability of Einstein spaces 
Sij = 0, can be seen quite easily if one decomposes the perturbation as 

hij = hf/' + Viipj + Vjcpi + f{x)gij, (69) 

which after integration by parts reduces the action (I55|) to 

^khXp - ( 70 ) 


where R is constant. The first derivative on the perturbation leads to a linear instability. This 
is consistent with our earlier flat space analysis. Hence Einstein spaces, including the flat space, 
critical points of the flow are linearly unstable. This is almost evident but let us see it more 
explicitly for the metric in the following coordinates 

ds"^ = cos^ p dr"^ + p dcj)^ + dp"^, (71) 

with R = 6. Assuming 


/ hii hi2 0 \ 

hX — hi2 —hii 0 , 

Vo 0 0/ 

second variation of the entropy about for this perturbation yields 

= _2 f Sshia, (72) 

Jm ^ 2 

which is not positive in general. Thus, is a saddle point not a minimum of the flow. In 0 , it 
was shown that as a critical point is stable under homogeneous deformations. 
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VI. EXTENDED COTTON SOLITONS 

In the previous section, we have shown that generic fixed points of Cotton flow are extended 
Cotton solitons that satisfy 


Cij + Vjlj' + VjVi + \{x)gij — 0, 


(73) 


where A(x) is not necessarily constant as was assumed in [13|. As was shown in 1^, if \(x) is 
constant and the manifold is compact without a boundary, then there are no non-trivial Cotton 
soliton solutions, namely A = 0, Cij = 0, and V) is a Killing vector 


£vg = 0 . 

For extended Cotton solitons, let us prove a similar theorem. 


(74) 


Theorem: On a compact Riemannian manifold without boundary, all solutions of (W 

have Cvg = —^g, namely, V is a conformal Killing vector and Cij = 0; There are no non-trivial 
extended Cotton solitons. 


Proof: The square of the Lie derivative of the metric reads 

\Cvg\^ = □(KiK*) - 2VinV^ + 2Vi{V^VjV^) - 2ViV^VjV^ - 2V^V^Rij. 

Trace and divergence of ([75]) can be computed as 


2ViK* + 3A(x) =0, 


(75) 


(76) 


□Ki + W^Vj - ^ViA(x) = 0. (77) 

Making use of these in d75|) . one arrives at 

\Cvg\^ = aiW) + 2Vi(K*A) + 2ViiV^VjV^) + 3X\x). (78) 

Integrating this over the manifold, once the boundary terms are dropped, yields 

J d^x^\Cvg\‘^ = ^ j d^Xy/g X‘^{x), (79) 

which can be written as 

J d^Xy/g\Cvg + X{x)g\'^ =0, (80) 

where we made use of (|70p . Since we are in a Riemannian manifold with a positive metric, from 
(j80|) it follows that V) is a conformal Killing vector with the conformal factor X{x). Then from (I73|) 
it follows that Cij = 0. So on a compact Riemannian manifold without a boundary, there are no 
non-trivial Cotton solitons besides the conformally flat ones. This proves the theorem. Of course 
for the Lorentzian signature and the non-compact manifolds, the theorem does not follow. 

Let us study some properties of extended Cotton solitons. Taking the divergence of ([77]) and 
making use of (I76p one arrives at 


-V'A(x)) =0, 


(81) 
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which can be recast as 


-d^X{x))) =0, (82) 

hence J* = — d^\) is a "conserved current". For some manifolds one can define the 

following conserved total charge. Snppose two space S foliate the three manifold Ai and let h* be 
the normal to the snrface S. And let jij denotes the indnced metric on the S. Then 

Q^J d'^x-y/yfii (^R^jV^ — d^X{x)J , (83) 

is a conserved total charge of the manifold, if the following condition is satisfied. Say we separate 
one of the coordinates as "r" and S as the constant r surface. Then for large r, one has 

R = R^.V^ -dA(x) - e>0, (84) 

large r r 


for i ^ r. This leads to 


^ = 0- 
dr 


(85) 


While, at the moment, we do know ex plic it examples of extended Cotton solitons, (one could 
perhaps find these nsing the method of [2^ since these are less restrictive solutions compared to 
the Cotton solitons with constant A). In fact let ns give an example in the Riemannian setting 
with a constant A 19(]. For constant A, there are many in the Lorentzian setting 19l42ll|. 

Consider the following metric 


( 1 \ ^ 
-{ydx — xdy) + dzj , 

with £ > 0. This metric has the following properties 


( 86 ) 


i) g = det[5rij] = 1 , 
iii)RijR^^ = 


a) R = — 




96 


(87) 


Also, the components of the Cotton tensor can be compnted as 


p £5^ ’ 

^ 482 

C22 = p-^^. 


Ci2 = -^xy, Ci^ = --y, 

^8 8 
C23 - -^x, C33 - 


( 88 ) 


This metric is a Cotton soliton for the scale factor A = —with the components of the vector 
field given as 


/ O ,0 iZ C2 \ 

~ + *^ 3 ^’ 'py + + '^22 -pZ - —X + Csy + caJ , 

where Cj’s are constants. For this metric, one can check that 

diiR^jV^) = 0, 


(89) 


(90) 


which is the A = constant version of (I82p . 








16 


It is quite interesting that (f55|) is also a Ricci soliton {Rij + ViXj + VjXi + Kgij = 0 ) with the 
scale factor A = -^ and the components of the vector field 

(X*) = (^ - - ciy + C 3 I, -^y + cix + C 2 , ^x + c^y + c^. (91) 

Using this fact we can define a new object that we shall call Topologically Massive solitons that 
solve 

Rij H —Cij + XiWj + VjWi + Rgij = 0, (92) 

l^ 

which is satisfied by (|86p with the following choices 

P = A + -A, Wi = Xi + -Vi. (93) 

/X /X 

Here g is the topological mass parameter. 


PP-Wave Solution of Cotton Flow 


In the Lorentzian setting, we can show that the pp-wave metric, (see 

ds'^ = 2dudv + dz^ + 2U{u,z)dv?, 


27| for its properties). 


(94) 


where u and v are null coordinates. The components of the metric can be written as = 
+ 2 UX^Xi, with the vector satisfying the following properties: 

A^A^ = 0, V^A, = 0, A^^S^U = 0, A^V,9^P = 0. (95) 

With these definitions, the Ricci tensor and the Cotton tensor can be computed as 

Pm- = -X^X.d^U, = e^\^X,)XaX^d^U. (96) 

The pp-wave metric (I94h satisfies the Cotton soliton equation with 

A = 0, (W) = (^0, CiXi + C3Z + C2, -C3XX - C4^ , 

^ pZ^y2(c3U+C4) -Zy/2{C3U+C4) (9^^) 

U(U, Z) = — ^ ' /l(tl) + ^ ' / 2 (») + flM, 

C3U + C4 y'2{c3U + C4) a/2(c3XX + C4) 

where Cj’s are arbitrary constants and fiiuYs are arbitrary functions. Interestingly, this Cotton 
soliton is also a gradient Cotton soliton [l9|] namely Cij + XiXjif: = 0 

^zb ^-zb 

i) = (p{u)+ lYz, U{u,z) =-^gi{u) - j^g 2 {u) + gz{u). (98) 

One can also show that the pp-wave metric is also a Ricci soliton with the following vector field 
Xi and the metric function U given as 

/ \ ca z e 

(X*) = (0, ciM + C32; + C2, -C3XX - C4), U{u,z) = -^-h-7-^- -hi{u)+ h2{u). ( 99 ) 

^ ' CsU + C4 2(C3Xi + C4) 

Finally, one can show that (|94D is a gradient Ricci soliton with 

^bz 

il; = ip(u)+ bz, U{u,z) =-^f{u) + g{u). 


( 100 ) 
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VII. CONCLUSIONS 

We have studied the linearized Cotton flow equation about a generic background with the help 
of a modified form of the DeTurck trick that removes the zero modes of the relevant operator in 
the flow equation. Then we specifically studied the perturbations about the critical points of the 
flow. We found that flat space and Einstein spaces, as critical points of the flow, have unstable 
modes at the perturbative level, making these spaces saddle points rather than minima. We have 
also shown that (in Appendix C) certain conformally flat fixed points that are not Einstein metrics, 
also have unstable modes with a dispersion relation cubic in the Fourier momentum. We have also 
supported our arguments by computing the second variation of the entropy functional about the 
critical points and show that the second variation is negative for Einstein metrics making them 
saddle points. Finding the true minima (if there is any) of the Cotton flow equations is an open 
problem. 

We have refined the gradient flow formulation of Cotton flow and gave a definition of extended 
Cotton solitons and worked out some properties of these critical points. We have also given an 
example of Topologically Massive solitons that is constructed from a solution to both Ricci and 
Cotton soliton equations. Our work could be relevant to both mathematics, as discussed in the 
Introduction and physics. For example it would be interesting to see the connection between the 
Cotton flow and various holographic theories that have the pure Chern-Simons action in the bulk 
of a three manifold [ 2 ^ . 


VIII. APPENDIX A: EVOLUTION OF SOME QUANTITIES TENSOR UNDER THE 

COTTON FLOW 


Under the Cotton flow, let us note how some tensors and scalars evolve, 

1. The Cotton tensor itself evolves as 

- +i^j. 


( 101 ) 


2. Square of the Cotton tensor evolves as 


+ 2CijRhV,C^m ) + 


3. The Riemann tensor evolves as 


^tR^,kl = + \{RkuPC,p - Rki/Cip) 


( 102 ) 


(103) 


4. Square of the Riemann tensor evolves as 


^tiR^JklR'^'^') = -^R'^aCij - lORRijC^^ + 16RjiR^^C\. 


(104) 
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5. Time independence of Bianchi identity dtiViC''^) = 0 leads to 


= 0 , = dtC^^ + - ^CkiO 


9'^^ ^ki 


(105) 


Note that at the critical points, the flow stops for all the above quantities. 


IX. APPENDIX B: LINEARIZATION OF THE COTTON TENSOR AROUND A 

GENERIC BACKGROUND 


Let us give some details of the calculations leading to equation Q. By perturbing the metric 
about an arbitrary background as gij = gij + hij, one obtains 

yk{Gh)L + rf'"' {TKn)L&i + i^3 + 0{h^). (106) 

A B 

From now on we will drop the 0{h?) terms. The background metric is compatible with the 
background covariant derivative ^kOij = 0, hence the linearized Christoffel connection reads 

(r-^fcn)L = -jg^'^iykhmn + ^nhmk “ ^mhkr^ ■ (107) 

Our task is to rewrite (|106l] in such a way that one can see the terms that can be related to 
diffeomorphisms of the background metric. For this purpose, let us first evaluate the A term in 
(llOOp . Note that the linearization of the Einstein tensor yields 

{G^i)l = Gml)L = -y^G^i + 9^^ {Gml)L. (108) 


Where {Gij)L = {Rij)L — ^gijRh — \hijR. By using the explicit forms of the Linearized Ricci 
tensor and Ricci scalar 


{Rij)L = - Uhij - 


Rl = {g^^Rij)L = -h^^Rij -ah + 
one can write the A part as 


(109) 


1 


Vfc(G^z)L = ( - ^kih^^Gml) + - ( VfcVnVffe"^. + - VkViV^h 


c 


D 


- ^VkihhR) 


Let us now calculate each term of (jllOp separately: the C term can be written as 

+ Vfc( - hRi^ + Wnh'^i - 

where we have used the three dimensional identity 

- ~ - R R 

Rkpqi ~ 9kqRip 9kiRqp T 9piRqk 9pqRik "^gkqgip “I 1^9ki9qp- 


( 110 ) 

( 111 ) 


(112) 
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Similarly the D and E terms can be written as 

v'’"' (v^VnVkhr - V^{R^khnl) - Rh^nhk" 

+ Vk(2KiR^^ - hRi^ + V^Rin - f ^'0) ’ 


( 113 ) 


VkDhh = {uVkhh + - RnkS^h^i) + {hJ^NnRK - h^N^Rnk) 

- Rl^Vnhh - 2R^k^lh^n + R^lhh - h^n^lR'k]- 


( 114 ) 


Using these results one arrives at 

7 ?*^' Vk{G^i)L = ( - 2hVkRi^ - RfVkh - h^iVkR + SR^^Vkhni + ^KNkR"^ 

+ Rnk^^hJ^l - h^NnRh + R^k^lh^n " Rkr^^hU ( 115 ) 

+ v^Vrykh^ -uVkV^. 


The B term in (jlOhp is easy to handle: 

{^Kn)LG^l = i?"z (Vfc V + - y^hun) ■ (116) 

Thus, collecting all these pieces, one arrives at the desired equation: 

2(C'^)l = - y UVkhh + V^VrS/kh^ + Vk{R^^Ki) - Ri^Vkh 

- h!^NnRh - VNkR + v'’"' Rnk^^h^i + R^^nkh + i^j. 

(117) 


X. APPENDIX C: LINEARIZED DEFORMATION OF THE CONFORMALLY FLAT 

BACKGROUND 

In this Appendix, let us give a more explicit linearization of the Cotton flow equation about its 
conformally flat fixed point. Any metric of the form 

= 0^(x, y, z) {dx^ + dy"^ + dz"^^ , (118) 

with a smooth 0 (x, y,z) > 0 is a critical point of the flow. Consider the most general perturbation 
about a given fixed point as 

ds^ = fl^(^(l + a) dx^ + 2bdxdy + 2cdxdz + (! + /) dy^ + 2gdy dz + {I + h) dz^'^, (119) 

where all the functions depend on x and t with the assumption that, save O, all functions vanish at 
large distances at t = 0. At first order in the perturbation theory, one computes the components 



20 


of the Cotton tensor as 

Cii = ^ (d^b - dydlc + d^dzh - dyC - d^d'^g - d^dyd^f + d^dyd^h + d^dyg^ , 

<^12 = ^ a + d^f - 2dydlg - d^dza + dyd^h + 2dxd‘^c + 2dxdyC + d^dj 

- dld^h - 2dldyg^ , 

Cis = ^ {dydla - dydlf + 2dydzg + d^a - dyh - 2dxd^b - 2dxdlb + 2dld^g 

+ dldyf - dldyh ^, ( 120 ) 

C 22 = ^ ( - Of 5 + dydlc + d^dlg + d^dyd^a - d^dyd^h - d^d^b - d^dyC + d^g^, 

C 23 = ^ i^dydlb - 2dydzC + d^d^a - d^dlf - d^d^a + d^d^h - 20 fO^c + 2dldyb 

- dxf + dxh ), 

C 33 = ^ ( - d^Ozb + dyC - dxdyd^a + d^dyd^f - d^d^g + Of0^6 + Of O^c - Of 5 ). 

Even though one can use the symmetries of the equation to reduce the system, this is still an 
unwieldy set. But our task is not to find the general solution rather to find the potentially unstable 
modes, hence let us assume for the sake of simplicity that the perturbations depend on 2 : and t but 
not on X and y. Then one arrives at 


dtin^a) = ^0f6, dtin^f) = -^0f6, 

dt{0,‘^c) = 0 , dt{VL^g) = 0 , 


dm^b) = ^dl{f-a), 

dt{n^h) = 0 . 


( 121 ) 


Further assuming c(0, z) = h(0, z) = g{0, z) = 0, the second line of (I12ip yields the vanishing of 
these functions for all t. From the first two equations of the first line, one finds f{t, z) = —a{t, z). 
Therefore, we have 

^ dtin'^b) = -^5fa, (122) 

which yield a linearized, complex KdV-type equation with a variable coefficient 

dt(n'^{a±ib)^ =^^d^(a±ib). (123) 

For flat space Q = 1 and the Fourier transform of this equation yields the modes found in the text 


u;(p) 



(124) 


with the + mode giving the perturbative instability. We have not been able to solve (| 12311 for 
generic Q{x,y, z,t), but assuming Q = i}{x,y) one finds the dispersion relation 


"■'P* = 

hence the unstable mode survives for this type of conformally flat background backgrounds. There¬ 
fore, it is clear that among the conformally flat fixed points of the flow, some metrics are saddle 
points rather than being the minima. It is an outstanding problem to find the minima of the 
Cotton flow. 
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XI. APPENDIX D: COTTON TENSOR UNDER AN ARBITRARY FLOW AND 

UNDER RICCI FLOW 


Our computation in Appendix B allows us to compute the flow of the Cotton tensor under an 
arbitrary geometric flow defined as 

dtQij = Eij, (126) 

where Eij is a symmetric tensor. With the help of equation (|117p . one gets 

2dtC^^ = _ M cij _ aVkE^i + Vk{R^^E^i) - R^VuE 

- E^iVr^Rh - EhVkR + RnkV^E^i + Ri^V^Eh + i^j. 

(127) 

It is interesting to see how the Cotton tensor behaves under the normalized and unnormalized 
Ricci flow which we give below. Cotton tensor under unnormalized Ricci flow dtQij = —2Rij flows 
as 

2dtC^^ = - 6R^^C\ + 2g^^RnkC'^^ + g^’^^Ri^VkR - 2rj^^^R^iVkR^^ +i^j- (128) 

Cotton tensor under normalized Ricci flow dtgij = —2{Rij — ^gijR) flows as 

28,c« = □£;« + (129) 

It is nice to see that both equations yield a parabolic flow of the Cotton tensor and hence amenable 
to maximum principle. 
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